The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom.

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Effects of relational and instrumental messaging on human perception of rattlesnakes.

We tested the effects of relational and instrumental message strategies on US residents' perception of rattlesnakes-animals that tend to generate feelings of fear, disgust, or hatred but are nevertheless key members of healthy ecosystems. We deployed an online survey to social media users (n = 1,182) to describe perceptions of rattlesnakes and assess the change after viewing a randomly selected relational or instrumental video message. An 8-item, pre-and post-Rattlesnake Perception Test (RPT) evaluated perception variables along emotional, knowledge, and behavioral gradients on a 5-point Likert scale; the eight responses were combined to produce an Aggregate Rattlesnake Perception (ARP) score for each participant. We found that people from Abrahamic religions (i.e., Christianity, Judaism, Islam) and those identifying as female were associated with low initial perceptions of rattlesnakes, whereas agnostics and individuals residing in the Midwest region and in rural residential areas had relatively favorable perceptions. Overall, both videos produced positive changes in rattlesnake perception, although the instrumental video message led to a greater increase in ARP than the relational message. The relational message was associated with significant increases in ARP only among females, agnostics, Baby Boomers (age 57-75), and Generation-Z (age 18-25 to exclude minors). The instrumental video message was associated with significant increases in ARP, and this result varied by religious group. ARP changed less in those reporting prior experience with a venomous snake bite (to them, a friend, or a pet) than in those with no such experience. Our data suggest that relational and instrumental message strategies can improve people's perceptions of unpopular and potentially dangerous wildlife, but their effectiveness may vary by gender, age, religious beliefs, and experience. These results can be used to hone and personalize communication strategies to improve perceptions of unpopular wildlife species.

Paradoxical Exception to Island Tameness: Increased Defensiveness in an Insular Population of Rattlesnakes.

Island tameness results largely from a lack of natural predators. Because some insular rattlesnake populations lack functional rattles, presumably the consequence of relaxed selection from reduced predation, we hypothesized that the Santa Catalina Island, California, USA, population of the southern Pacific rattlesnake (Crotalus helleri, which possesses a functional rattle), would exhibit a decrement in defensive behavior relative to their mainland counterparts. Contrary to our prediction, rattlesnakes from the island not only lacked tameness compared to mainland snakes, but instead exhibited measurably greater levels of defensiveness. Island snakes attempted to bite 4.7 times more frequently as we endeavored to secure them by hand, and required 2.1-fold more time to be pinned and captured. When induced to bite a beaker after being grasped, the island snakes also delivered 2.1-fold greater quantities of venom when controlling for body size. The additional venom resulted from 2.1-fold larger pulses of venom ejected from the fangs. We found no effects of duration in captivity (2-36 months), which suggests an absence of long-term habituation of antipredator behaviors. Breeding bird surveys and Christmas bird counts indicated reduced population densities of avian predators on Catalina compared to the mainland. However, historical estimates confirmed that populations of foxes and introduced mammalian predators (cats and pigs) and antagonists (herbivorous ungulates) substantially exceeded those on the mainland in recent centuries, and therefore best explain the paradoxically exaggerated defensive behaviors exhibited by Catalina's rattlesnakes. These findings augment our understanding of anthropogenic effects on the behaviors of island animals and underscore how these effects can negatively affect human safety.

Distinct cardiotoxic effects by venoms of a spitting cobra (Naja pallida) and a rattlesnake (Crotalus atrox) revealed using an ex vivo Langendorff heart model.

Here we describe the acute myocardial effects of an elapid (red spitting cobra, Naja pallida) and a viper (western diamondback rattlesnake, Crotalus atrox) venom using an ex vivo heart model. Our results reveal two different pathophysiological trajectories that influence heart function and morphology. While cobra venom causes a drop in contractile force, rattlesnake venom causes enhanced contractility and frequency that coincides with differences in myocellular morphology. This highlights the medical complexity of snake venom-induced cardiotoxicity.

Comparative analysis of Deinagkistrodon acutus venom from Taiwan and China utilizing chromatographic, electrophoretic, and bioinformatic approaches, along with ELISA employing a monospecific antivenom.

Deinagkistrodon acutus is a medically important pitviper inhabiting mainly South China and Taiwan. The hemorrhagic effects of its envenoming are compatible to its venom, which is abundant in metalloproteases (svMPs) and C-type lectin-like proteins. In this study, we investigated geographic variations in the venom of D. acutus collected from Taiwan and four Mainland Chinese provinces: Fujian, Jiangxi, Anhui, and Hunan. The variations were assessed through high-performance liquid chromatography, non-metric multidimensional scaling analysis, gel electrophoresis, and enzyme-linked immunosorbent assay (ELISA) with a monospecific antivenom (DaMAV) generated against the Taiwanese D. acutus venom, and discussed based on venom-protein sequences in databases and literature related to D. acutus venom. Additionally, the cross-reactivity of DaMAV against Crotalus horridus and Calloselasma rhodostoma venoms was investigated. We noted differential abundances of D. acutus venom metalloproteases, C-type lectin-like proteins, and phospholipase A2, along with point mutations and selective expression of serine protease isoforms. The ELISA results revealed that the venom from Taiwan was more reactive toward Taiwanese DaMAV than the four Mainland Chinese venoms, consistent with chromatographic profile differences, whereas C. horridus venom presented moderate cross-reactivity with DaMAV. The observed immunoreactivities of these venom with DaMAV can be attributed to the high prevalence of their PIII-svMPs, which are the dominant antigens, and the conservation of PIII-svMP epitopes.

Motor control: Snake neurons speed up.

How are motor neurons tuned for very different jobs? Classic work has focused on variations in motor neuron size and their premotor networks. New results in rattlesnakes show that shifting a motor neuron's temporal precision can be as simple as changing its potassium channel conductance.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes.

BACKGROUND: Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris). RESULTS: Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis. CONCLUSIONS: Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Effects of the Heterodimeric Neurotoxic Phospholipase A2 from the Venom of Vipera nikolskii on the Contractility of Rat Papillary Muscles and Thoracic Aortas.

Phospholipases A2 (PLA2s) are a large family of snake toxins manifesting diverse biological effects, which are not always related to phospholipolytic activity. Snake venom PLA2s (svPLA2s) are extracellular proteins with a molecular mass of 13-14 kDa. They are present in venoms in the form of monomers, dimers, and larger oligomers. The cardiovascular system is one of the multiple svPLA2 targets in prey organisms. The results obtained previously on the cardiovascular effects of monomeric svPLA2s were inconsistent, while the data on the dimeric svPLA2 crotoxin from the rattlesnake Crotalus durissus terrificus showed that it significantly reduced the contractile force of guinea pig hearts. Here, we studied the effects of the heterodimeric svPLA2 HDP-1 from the viper Vipera nikolskii on papillary muscle (PM) contractility and the tension of the aortic rings (ARs). HDP-1 is structurally different from crotoxin, and over a wide range of concentrations, it produced a long-term, stable, positive inotropic effect in PMs, which did not turn into contractures at the concentrations studied. This also distinguishes HDP-1 from the monomeric svPLA2s, which at high concentrations inhibited cardiac function. HDP-1, when acting on ARs preconstricted with 10 µM phenylephrine, induced a vasorelaxant effect, similar to some other svPLA2s. These are the first indications of the cardiac and vascular effects of true vipers' heterodimeric svPLA2s.

COMPARING THE EFFECTS OF LITHIUM HEPARIN AND DIPOTASSIUM ETHYLENEDIAMINETETRAACETIC ACID ON HEMATOLOGIC VALUES IN PRAIRIE RATTLESNAKES (CROTALUS VIRIDIS) AND LAKE ERIE WATER SNAKES (NERODIA SIPEDON INSULARUM).

Anticoagulants prevent clotting of blood samples and preserve cellular morphology for hematologic evaluations, but studies comparing anticoagulants are limited in snakes. The objective of this study was to evaluate the effects of lithium heparin (LH) and dipotassium ethylenediaminetetraacetic acid (EDTA) on hematologic values in prairie rattlesnakes (PR; Crotalus viridis, n = 16) and Lake Erie water snakes (LEWS; Nerodia sipedon insularum, n = 21). Venipuncture was performed and blood samples were immediately aliquoted into LH and EDTA microtainers. Packed cell volume (PCV), total solids (TS), 100-cell differential counts, and Avian Leukopet white blood cell counts (WBC) were determined for each anticoagulant. Passing-Bablok regression and Bland-Altman plots revealed that anticoagulant choice did not constantly or proportionally bias the values of any WBC parameter. Mixed models demonstrated that blood anticoagulated with EDTA had higher PCV in PR (P = 0.04) and TS in both species (P < 0.05). However, the magnitude of the differences attributable to anticoagulant choice was relatively small and likely not clinically important. Hemolysis was not appreciated in any samples. Our findings demonstrate that LH and EDTA are equally appropriate for use in PR and LEWS, but may require separate reference values.

Venom of the neotropical rattlesnake, Crotalus culminatus: Intraspecific variation, neutralization by antivenoms, and immunogenicity in rabbits.

Crotalus culminatus is a medically significant species of rattlesnake in Mexico [1]. While the proteomic composition of its venom has been previously reported for both juvenile and adult specimens, there has been limited research into its functional properties, with only a few studies, including one focusing on coagulotoxicity mechanisms. In this study, we aimed to compare the biochemical and biological activities of the venom of juvenile and adult snakes. Additionally, we assessed antibody production using the venoms of juveniles and adults as immunogens in rabbits. Our findings reveal lethality and proteolytic activity differences between the venoms of juveniles and adults. Notably, juvenile venoms exhibited high proportions of crotamine, while adult venoms displayed a reduction of this component. A commercially available antivenom demonstrated effective neutralization of lethality of both juvenile and adult venoms in mice. However, it failed to neutralize the paralytic activity induced by crotamine, which, in contrast, was successfully inhibited by antibodies obtained from hyperimmunized rabbits. These results suggest the potential inclusion of C. culminatus venom from juveniles in commercial antivenom immunization schemes to generate antibodies targeting this small myotoxin.

Genetic mechanisms and biological processes underlying host response to ophidiomycosis (snake fungal disease) inferred from tissue-specific transcriptome analyses.

Emerging infectious diseases in wildlife species caused by pathogenic fungi are of growing concern, yet crucial knowledge gaps remain for diseases with potentially large impacts. For example, there is detailed knowledge about host pathology and mechanisms underlying response for chytridiomycosis in amphibians and white-nose syndrome in bats, but such information is lacking for other more recently described fungal infections. One such disease is ophidiomycosis, caused by the fungus Ophidiomyces ophidiicola, which has been identified in many species of snakes, yet the biological mechanisms and molecular changes occurring during infection are unknown. To gain this information, we performed a controlled experimental infection in captive Prairie rattlesnakes (Crotalus viridis) with O. ophidiicola at two different temperatures: 20 and 26°C. We then compared liver, kidney, and skin transcriptomes to assess tissue-specific genetic responses to O. ophidiicola infection. Given previous histopathological studies and the fact that snakes are ectotherms, we expected highest fungal activity on skin and a significant impact of temperature on host response. Although we found fungal activity to be localized on skin, most of the differential gene expression occurred in internal tissues. Infected snakes at the lower temperature had the highest host mortality whereas two-thirds of the infected snakes at the higher temperature survived. Our results suggest that ophidiomycosis is likely a systemic disease with long-term effects on host response. Our analysis also identified candidate protein coding genes that are potentially involved in host response, providing genetic tools for studies of host response to ophidiomycosis in natural populations.

Structural studies with crotoxin B from Crotalus durissus collilineatus venom suggest a heterodimeric assembly formed by two new isoforms.

In accidents involving Crotalus snakes, the crotoxin complex (CTX) plays lethal action due to its neurotoxic activity. On the other hand, CTX have potential biotechnological application due to its anti-tumoral, anti-inflammatory, antimicrobial, analgesic and immunomodulatory properties. CTX is a heterodimer composed of Crotoxin A (CA or crotapotin), the acidic nontoxic and non-enzymatic component and; Crotoxin B (CB), a basic, toxic and catalytic PLA2. Currently, there are two classes of CTX isoforms, whose differences in their biological activities have been attributed to features presented in CB isoforms. Here, we present the crystal structure of CB isolated from the Crotalus durissus collilineatus venom. It amino acid sequence was assigned using the SEQUENCE SLIDER software, which revealed that the crystal structure is a heterodimer composed of two new CB isoforms (colCB-A and colCB-B). Bioinformatic and biophysical analyses showed that the toxin forms a tetrameric assembly in solution similar to CB from Crotalus durissus terrificus venom, despite some differences observed at the dimeric interface. By the previously proposed classification, the colCB-B presents features of the class I isoforms while colCB-A cannot be classified into classes I and II based on its amino acid sequence. Due to similar features observed for other CB isoforms found in the NCBI database and the results obtained for colCB-A, we suggest that there are more than two classes of CTX and CB isoforms in crotalic venoms.

Unveiling the secrets of snakes: Analysis of environmental, socioeconomic, and spatial factors associated with snakebite risk in Paraná, Southern Brazil.

The state of Paraná is home to three out of the five medically significant snake genera in Brazil and lacks of snakebite epidemiology studies. This study aimed to ascertain the spatial, environmental, and socioeconomic factors associated with snakebite risk by analyzing notification data of cases in the state of Paraná. Notification and socioeconomic data were gathered from the online platforms of the National System of Notifiable Diseases (SINAN) and the Brazilian Institute of Geography and Statistics (IBGE). Land cover and land use maps were obtained from the Mapbiomas platform in raster format and subsequently converted into vectors using QGis software. The proportions of land use and land cover in square kilometers (km2) were then calculated. All acquired data were tabulated using Microsoft Excel 365 software. For spatial analysis, GeoDa software version 1.20 was utilized to calculate the Global and Local Moran indices, assessing spatial correlations. Between 2007 and 2021, 12,877 notifications were recorded, with an average incidence of 8.22/100,000 inhabitants in the state, 8166 (63.41%) caused by Bothrops, 1534 (11.91%) caused by Crotalus, 56 (0.43%) caused by Micrurus. 1703 (13.22%) caused by non-venomous snake species, and the remaining cases did not have the identified causative species. The incidents caused by Bothrops and Crotalus showed different distribution patterns. Spatial analysis revealed that key factors contributing to snakebite risk included the presence of native forests, mangroves, apicuns, and monospecific planted forests. The population group at the highest risk comprised rural residents and workers. Furthermore, the absence of basic sanitation and proper garbage collection and disposal exhibited positive correlations with snakebites. Conversely, intensive farming practices with substantial mechanization and pastures demonstrated negative spatial correlations. This study has enabled the identification of the primary factors associated with snakebite risk, facilitating more targeted efforts to prevent snakebite accidents among vulnerable populations.

Timing and precision of rattlesnake spinal motoneurons are determined by the KV72/3 potassium channel.

The evolution of novel motor behaviors requires modifications in the central pattern generators (CPGs) controlling muscle activity. How such changes gradually lead to novel behaviors remains enigmatic due to the long time course of evolution. Rattlesnakes provide a unique opportunity to investigate how a locomotor CPG was evolutionarily modified to generate a novel behavior-in this case, acoustic signaling. We show that motoneurons (MNs) in the body and tail spinal cord of rattlesnakes possess fundamentally different physiological characteristics, which allow MNs in the tail to integrate and transmit CPG output for controlling superfast muscles with high temporal precision. Using patch-clamp electrophysiology, we demonstrate that these differences in locomotor and rattle MNs are mainly determined by KV72/3 potassium channels. However, although KV72/3 exerted a significantly different influence on locomotor and rattle MN physiology, single-cell RNA-seq unexpectedly did not reveal any differences in KV72/3 channels' expression. VIDEO ABSTRACT.

Locomotor kinematics on sand versus vinyl flooring in the sidewinder rattlesnake Crotalus cerastes.

For terrestrial locomotion of animals and machines, physical characteristics of the substrate can strongly impact kinematics and performance. Snakes are an especially interesting system for studying substrate effects because their gait depends more on the environment than on their speed. We tested sidewinder rattlesnakes (Crotalus cerastes) on two surfaces: sand collected from their natural environment and vinyl tile flooring, an artificial surface often used to elicit sidewinding in laboratory settings. Of ten kinematic variables examined, two differed significantly between the substrates: the body's waveform had an average of ∼17% longer wavelength on vinyl flooring (measured in body lengths), and snakes lifted their bodies an average of ∼40% higher on sand (measured in body lengths). Sidewinding may also differ among substrates in ways we did not measure (e.g. ground reaction forces and energetics), leaving open clear directions for future study.

Antileishmanial effects of γCdcPLI, a phospholipase A2 inhibitor from Crotalus durissus collilineatus snake serum, on Leishmania (Leishmania) amazonensis.

BACKGROUND: Leishmaniasis, a neglected disease caused by the parasite Leishmania, is treated with drugs associated with high toxicity and limited efficacy, in addition to constant reports of the emergence of resistant parasites. In this context, snake serums emerge as good candidates since they are natural sources with the potential to yield novel drugs. OBJECTIVES: We aimed to show the antileishmanial effects of γCdcPLI, a phospholipase A2 inhibitor from Crotalus durissus collilineatus snake serum, against Leishmania (Leishmania) amazonensis. METHODS: Promastigotes forms were exposed to γCdcPLI, and we assessed the parasite viability and cell cycle, as well as invasion and proliferation assays. FINDINGS: Despite the low cytotoxicity effect on macrophages, our data indicate that γCdcPLI has a direct effect on parasites promoting an arrest in the G1 phase and reduction in the G2/M phase at the highest dose tested. Moreover, this PLA2 inhibitor reduced the parasite infectivity when promastigotes were pre-treated. Also, we demonstrated that the γCdcPLI treatment modulated the host cell environment impairing early and late steps of the parasitism. MAIN CONCLUSIONS: γCdcPLI is an interesting tool for the discovery of new essential targets on the parasite, as well as an alternative compound to improve the effectiveness of the leishmaniasis treatment.

Early-Onset Respiratory Muscle Paralysis in Crotalic Envenomation: A Case Study.

Crotalic envenomation is responsible for approximately 8%-13% of ophidism cases in Brazil, yet it is associated with the highest mortality among snakes. We describe the case of a patient bitten by a rattlesnake who developed ventilatory muscle paralysis within hours after envenomation. While diaphragmatic paralysis is a rare late neurotoxic event following crotalic envenomation, in this case, paralysis occurred early but was rapidly reversed after antivenom administration. This report discusses potential contributing factors based on a comprehensive literature review.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake.

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.